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X ■ 1. INTRODUCTION 

^' 

Binary star systems are frequently the only source for the fundamental determination of the 
most basic property of a star, its mass. The masses of a sample of stars, when determined with 
sufficient accuracy, serve as crucial tests of our theoretical understanding of stellar formation, 
structure, and evolution. By combining information from the angular orbit projected on the sky 
(or from a light curve in the case of an eclipsing system) with the spectroscopic orbit, one can 
determine the distance and masses of the individual stars. In the cases of the most massive stars, 
accurate masses are even more important, not only because massive stars are relatively rare, but 
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because these stars (initial masses > 20 M(7)) play critical ro les in galaxies. Massive stars serve 
as signposts of star formation in galaxies (jMassev et al.lll995l ). Current theory suggests that high 
mass stars are born in loosely bound star clusters which disperse after a few Myr. These birth 
conditions su ggest that high ma ss stars might often have companions, as for example, found around 
7^ Velorum ( Pozzo et al.ll200d ). Some of these could be very close and gravitationally bound to 
the O star, which is suggested by the prevalence of binary and triple systems among young O stars 
(|Mason et alJIlQQsl V 



Mason et al.l (|l998l ) performed a speckle interferometry survey of Galactic 0-type stars for close 



companions, specifically looking for differences in the multiplicity frequencies amongst the cluster, 
field, and runaway 0-type star populations, as well as the distribution of orbital periods. They 
did their survey in the T^-band using speckle interferometry, which is sensitive to the detection of a 
companion if the projected separation is in the range 0'.'035 < p < l''5 and the magnitude difference 
Amy < 3. Our survey is in the I-band (which will slightly emphasize redder companions) and 
extends the dynamic range of their work. Due to the availability of only one photometric band 
for this survey, we are unable to determine physical associations for the new companions with 
any certainty, and we can only address the likelihood that perhaps some of the companions are 
gravitationally bound. The results presented here invite further investigation into the O stars and 
their potential to shed light on star formation processes. 



2. OBSERVATIONS 

The data were taken using the adaptive optics (AO) system and Visible Imager (Visim) camera 
of the Advanced Electro-Optical System (AEOS) 3.6-m telescope at the Maui Space Surveillance 
System, located on Haleakala. The data were collected during four separate observing runs between 
2001 February and 2002 September. In addition to the dedicated observing runs, some observa- 
tions were taken as part of a queue-scheduled observation program between 2001 May and 2005 
November. 



2.1. AEOS Telescope and Adaptive Optics System 



The AEOS telescope is an altitude-azimuth (alt- az) design. While the teles cope is very flexible 
and supports several different optical configurations (iRoberts &: Neymanll2002l ). of interest here is 
the Cassegrain configuration which feeds a 726-m focal length, //200 beam (with a field of view 
of 61'.'9) through a Nasmyth port coude path to a fold flat below the telescope which directs the 
beam onto a custom optical table that supports the main components of the AO system. For 
this project, the short-wavelength light (400 ~ 700 nm) is sent to the wavefront sensor (WFS) and 
tracker camera beam trains, while the 700 nm to 5 pm light is sent to the VisIm, which is used as 
the science detector. 
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The heart of the AO system is a Shack-Hartmann WFS driving a 941-actuator Xinetics de- 
formable mirror (DM) with a stroke of ±4//m. Optics inside the WFS produce a pupil image on a 
Hartmann lens array. Each generated Hartmann spot is then imaged onto a 4x4 group of pixels 
on the WFS CCD. These pixel values can be used directly or binned (for higher throughput) to 
generate the wavefront slopes. The WFS CCD is a Lincoln Labs frame-transfer device with 16 
output ports, capable of frame rates of 0.2 to 2.5 kHz. The wavefront slopes are calculated from 
the Hartmann spot values by means of many digital signal processors working in parallel, feeding 
these data to a system which reconstructs the wavefront, taking into account alt-az induced image 
rotation and WFS-to-DM-actuator mapping artifacts. 



The Visim camera is described in detail in [Roberts &: NevmanI (|2002l ). In short, it operates 



from 700 to 1050 nm; is equipped with an atmospheric dispersion corrector; has a two-mode image 
derotator (zenith at a fixed position in the image or celestial north at a fixed position in the image); 
and, for this project, has a 10" field of view (0'.'020 pixel"^). The detector is a 512x512 pixel EEV 
CCD with a dark current of 22 e~ pixel"^ s~^ (at —40° C) and a read noise of 12 e" rms. The 
camera output is digitized to 12 bits with 10 e~ per digital number. 



2.2. Object Selection 



The object list started out as all the components in the 228 systems listed in iMason et al 



(jl998l ). They were then culled for the effective magnitude limit of the AEOS AO system (about 
my = 8), and for the declination limit of acceptable AO correction (objects that at some time 
during the year get above 30° elevation at Haleakala, i.e., a declination greater than about —45°). 
This red uced the list to 164 object s. We added seven more targets from the Galactic O Star 
Catalog ( Maiz-Apellaniz et al. 2004) that are within our adopted magnitude and declination limits. 
We actually observed, at least once, 116 of the 171 objects in our list. These observations are 
fairly consistently distributed among the three 0-star populations, sampling 63% of the cluster 
membership stars, 78% of the field stars, and 88% of the runaway stars in our list of potential 
targets. 



2.3. Data Collection and Reduction 

Since VisIm is only a 12-bit camera, a stellar image can overflow the digitizer in a rather 
short exposure time. If the frame saturates, then it is omitted from the final image. Frequently, 
this saturated frame will represent a period of particularly good seeing. In order to keep as many 
of these "good seeing" frames as possible, we set the exposure time such that the average peak 
value was about 75% of the full-well depth. We then built up the signal-to-noise ratio by taking 



^http: //www-int . stsci . edu/^jmaiz/reseajch/GDS/GOSmain.html 
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many frames, then weighting and summing them. Figure [T] shows the effect of this summing on 
the detectable magnitude difference for 1, 10, 50, 100, 250, and 1000 frame(s). On the basis of this 
plot, we aimed to take at least 1000 frames of each object, although time constraints, weather, and 
object brightness occasionally limited us to a lower frame count. 

The Visim is a frame transfer camera and is therefore unable to take a bias frame. To com- 
pensate, a bias frame was created by taking many dark exposures at a variety of exposure times 
ranging from 10 ms to several minutes. A linear fit was made to the signal from each pixel as a 
function of exposure time; the y-intercepts became the bias frame while the slopes became the dark 
frame. Prior to 2005, fiat - field f rames were created using the twilight sky to evenly illuminate the 
VisIm detector. iRoberts ( 2001 ) found that sky flat-field frames varied at about the 1% level, so 
after 2005, flat-field frames were generated using an internal calibration sphere; this lowered the 
frame-to-frame variation below the 0.01% level. From this point, all data frames were debiased, 
dark-subtracted and flat-fielded in the conventional way. 

For a given sequence of data frames, a weighted shift-and-add algorithm was used to create the 
final image. The weighted shift-and-add algorithm (jten Brummelaar et al.lll998l ) is a modification 
of the traditional image-stacking algorithm that takes seeing conditions in each individual frame 
into account. Frames with higher peak values (which represent better seeing) influence the final 
image (by means of weighting) more than frames with lower peak values. 



The fitting algorithm is described in detail in Iten Brummelaar et al.l (|l996l . |2000|). In short, 
the point-spread function (PSF) used to represent the system performance is that of the primary 
star in the image, with a few modifications; near the primary, the PSF is a pixel-for-pixel table of 
numbers, while farther out the PSF is a radially-symmetric series of values. This PSF is fitted to 
the primary and secondary components and iterated until the intensity ratio converges. 



2.4. Detection Sensitivity 

Figure [1] only gives a qualitative sense of the companion detection sensitivity of the system. 
Cross comparisons of the reduced images, especially of the same object taken on different nights, 
show different detection limits. This is primarily due to variations in AO performance which is 
primarily attributed to variations in the atmospheric turbulence profile. In order to quantify the 
sensitivity o f the reduced irnages , we have created a variation of the "dynamic range map" technique 



described in lHinklev et al.l (|2007l ) which defines the dynamic range of a given position in a 2-D image 
as the faintest companion detectable at that position to the 5a level. In our version, we construct 
a map the same size as each reduced image. The intensity level of each pixel in the map is set to 
five times the RMS intensity variation across a patch centered on the corresponding pixel in the 
original image. The patch is a square with lengths equal to the FWHM of the original image. This 
produces the dynamic range map in intensity terms, which are then converted to magnitudes. 

Figure [2] shows an example of this technique applied to one of the images generated in this 
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Radius (arcsec) 



Fig. 1. — Assuming that a Gaussian-shaped peak, 3a detection above background constitutes a 
positive detection of a companion, this plot shows the rough dynamic range of the AEOS AO 
system as a function of radius for different numbers of summed images. Prom top to bottom, the 
curves represent the results of 1000, 250, 100, 50, 10, and 1 summed frame(s). The "shelf", most 
apparent in the upper curves from about 0.7 to 1 arcsec in radius, is due to AO system artifacts. 
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survey. It is apparent from the sample figure that the detection threshold for a companion is highly 
spatially variant. The dynamic range increases with increasing radius from the central star. There 
are several artifacts in this image which lower the dynamic range. The largest is the vertical line, 
which is an artifact of the shutter-less frame tran sfer process. There ar e also four waffle spots caused 



by the Fried geometry of the wavefront sensor (jMakidon et al.ll2005l ). There are also hints of the 



diffraction pattern caused by the secondary mirror support spiders. As the AO system performance 
decreases, the FWHM of the PSF will increase, which increases the area over which the central PSF 
causes confusion. In addition, it also widens the companion's PSF, lowering the contrast between 
the two and decreasing the detection rate. After 2003, the waffle spots were greatly reduced when 
new reconstructors were put in place to filter out the unsensed waffle modes. 

These maps have several purposes. For those objects with multiple observations where com- 
panions are seen in some but not all of the data, changes in the dynamic range can help explain 
the reasons for these discrepancies. They also illustrate the limitations of the data, showing where 
it is almost impossible to find faint companions. This is often useful when trying to detect known 
long-period spectroscopic companions. 



3. RESULTS 

3.1. New Pairs and New Measures 

Table [1] lists the 40 new and 24 known pairs detected in the survey. Four of the new compan- 
ions (HD 17505, HD 37468, HD 193793, and HD 217086) are in systems with previously known 
components in the AEOS field of view. The remainder are new companions in systems not pre- 
viously known to be multiple (at least in the AEOS field of view). The table has nine columns 
which summarize the information gleaned from the measurements: Column (1) lists the Washing- 
ton Double Star (WDS) identifier, columns (2) and (3) identify the object by HD and Hipparcos 
number, and column (4) gives the discoverer designation (TRN for new pairs). Circumstances of 
the observation are given in columns 5-8. Column (5) gives the Besselian epoch of the start of the 
observation, columns (6) and (7) give the position angle (in degrees, where north is 0° and east 
is 90°) and separation (in seconds of arc) between the primary (the survey star - the brightest 
star in the image) and the detected component. Notice that for many of the objects, the position 
angle is missing or ambiguous. The Visim derotator is usually in one of two modes, one which 
keeps zenith up in the image, and one which keeps north up in the image. In the observations with 
more than one or a missing position angle, the state of the derotator was unknown, leaving the 
orientation of the image on the sky unknown. However, in some cases, position of a known pair 
resolves the ambiguity while for the rest, new observations will be required. Column (8) gives the 
magnitude difference at I-band between the primary and the detected component. Finally, column 
(9) lists the spectroscopic status. Given the separations detected here, none of the short period 
spectroscopic pairs are expected to be detected. However, listing the spectroscopic binaries helps 
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Fig. 2. — A sample dynamic range map (in this case, HD 34656, taken Bess. Year 2003.7360). The 
gradient map on the right shows the 5a detection hmit for a given pixel of the image in magnitudes. 
There are several artifacts in this image which lower the dynamic range. The largest is the vertical 
line, which is an artifact of the shutter-less frame transfer process. There are also four wafHe 
spots caused by the Fried geometry of the wavefront sensor. There are also hints of the diffraction 
pattern caused by the secondary mirror support spiders. These all increase the difficulty of finding 
any companions that coincide with these artifacts. 
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to form the complete multiplicity picture. The indicator in this table is V (velocity variable) if 
one or more close spectroscopic companions is indicated or C (velocity constant) if r adial velocity 



work indicates no close companion. The spectroscopic reference is generally given in lMason et al, 
(|l998l l. 



The differential magnitude errors in column 8 of Table [U were generated from simulated binary 
data that we analyzed using the same fitting algorithm applied to the measured binary data. 
To generate the simulated binaries, we started by collecting all the single star images from this 
survey (see Table H]) and those of a duplicity survey B stars carried out contemporaneously (see 
Roberts et al.l 120071 . Table 2). We omitted single star images where there were fewer than 250 
frames used to make the final image. This exercise gave us 167 unique images with PSFs ranging 
from 0'.'083 to 0'.'445, with a decided emphasis on values less than 0'.'2. For each entry in Tabled! 
167 simulated binaries were constructed using the magnitude difference in column 8 and the pixel 
positions of the primary and secondary components in the main image. We then ran the fitting 
algorithm on each of these 167 simulated images for each table entry. We calculated the magnitude 
variance in two passes, the first to filter out outliers on the basis of the variance of the fitted pixel 
separation of the simulated binary (namely, any deviation greater than la), and the second to use 
this reduced list of simulated fits to generate the final magnitude variance. In general, the number 
of outliers per table entry was less than 10. The errors quoted in column 8 are the standard errors of 
these variances. Errors are not quoted for the separation (column 7), but in all cases, the standard 
error is less thar i 0^.^04, usually less th an O'.'OOS. For errors in position angle (column 6), we adopt 
the values from [Roberts et al.l (|2007l ). ±2° for separations less than 1", and ±1° for separations 
greater than l" . Two effects lead to this uncertainty, the clocking calibration of the dove prism in 
the derotator, and the slight natural image rotation during the frame collection sequence when the 
derotator is in the fixed zenith position mode. 



3.1.1. Notes to Pairs 

Here we discuss the various binary systems in Table [J by WDS identifier and discovery desig- 
nation. In the cases of the stars classified as runaways, it is unlikely that the detected companion 
is physical (see § H]) . 

02229+4129 = TR N 10: This is a runaway star with a close, low-luminosity companion 
(jBoyajian et al.ll2005l ). 

02407+6117 = 

0251 1+6025 = TRN 13 AH: The A component of HD 17505 is a spectroscopic triple (iHillwig et al, 
20061 ). It is uncertain if the new component (AH, p ~ 4'.'6) and the known pair of Table [1] 



TRN 12AD: The A component is a spectroscopic triple (jMcSwainI l2003l ) . 



(STF 306AB ) are both gravitationally bound. If STF 306AB is indeed gravitationally bound, 
Hillwig et al.l (|2006l ) estimate it to have a period of 27,000 years. 
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03556+5238 = HDS 494: The Am was probably a little too la rge for iMason et al.1 (Il998ll to 



detect this pair, but it was detected in later speckle observations by IMason et alj (|2001a 
03590+3548 = TRN 16: Classified as a runaway star in Icies BoltonI (|l98d ). 



05163+3419 = TRN 17Aa: Classified as a runaway star in lCies BoltonI (11989 ). SET 136AB 
(p ~ S'.'S, Am ~ 3.3, just outside the field of view) is a likely optical pair, with the B component 
not physically associated with this close double. 

05207+3726 = TRN 18: One new close pair and three wide pairs are found here. The close one 
is designated Aa while the wider pairs are designated AB, AC, and AD. The known pair SEX 201 
(p 24'.'7, Am ~ 5.9) is outside the field of view and has been designated AE. 

05297+3523 = HU 217: This known, relatively close, pair has been measured since 1900. The 
position angle has only changed by 4° over this time span. 

05353-0523 = STF 748Ca,F: The primary is 9^ Ori C, the mos t massive star in the Orion 
Trapezium. It has a closer companion discovered by speckle methods (jKraus et al.l 120071 ) . 



05354-0525 = CHR 24 9Aa: HIPPARCOS (IPerrvman fc ESAlll997l ) measured a AHp of 3.23. 
Fabricius fc Makarovl (|200nl ) found ABt = 4.37 ± 0.02 and AVr = 3.21 ± 0.01. 



05387—0236 = TRN 19AF: The new component (AF) is at a separation nicely nestled between 
BU 1032AB and STF 762AC, at a value which would be consistent with a hierarchical arrangement. 
The new TRN 19AF compo nent is much brighter in the mid-infrared where it was first discovered 



van Loon Oliveiral l2003l ) and it is probably a low mass K-star surrounded by a circumstellar 



disk (jCaballerQ .20071 ) . See § 13.1.31 below for a discussion of the orbit of BU 1032AB. 

05407-0157 = STF 774Aa-B: This pair has an orbit (|HoDmanJl"967l ') with residuals of: 0-C = 
—5° and -^.'03 in position angle and separation. However, the orbit is classified as "indeterminate" 
(grade 5) in the 6th Catalog of Orbits of V.suai B.nary StaS , SO these residuals are not indicative 
of measurement quality. While the total number of measures has increased by about 50% since this 
orbit was published, the initial elements are of such low quality that even at this point any sort of 
improvement above "indeterminate" subjective quality is not yet possible. Some of the scatter in 
the residuals may be due to the presence of a closer companion to C Ori A found by interferometry 
(|Hummel et al.lboool ). 



06319+0457 = GAN 3AB: The m easure agrees with earlier published data quite well and 



indicates that the measure made in 1928 (jSilbernagellll93ll ) is probably erroneous. 



06410+0954 = STF 950 Aa-B: This pair has only changed by ~9° and 0'.'2 since it was first 
resolved by Struve in 1825 ()Struvdll837l ). 



16550-4109 = I 576: I 576 has not been measured since 1934 (jWallenquist] 1 19341 ). so 



^http : //ad . usno . navy . mil/wds/orb6 . html | 
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identification of this pair with that of Innes is uncertain. 



17065-3527 



B 



894: This binary, unconfirmed since its first resolution by Ivan den Bos 



(119281 ) ■ has closed slightly. 

17158—3344 = SEE 322: This pair h as showri noticeable Keplerian motion, completing about 
70° of orbital revolution since discovery (jSedll898l ). See S 13.1.31 below for a discussion of the orbit. 

17347-3235 = HDS2480Ab: The pairs, ISO 5 and HJ 4962 should have been seen as well. 
ISO 5 was probably missed due to the corrected FWHM. HJ 4962 was probably just out of the 
field. 

18024-2302 = TRN 26AH: While this pair may be physical, the known double H N 40 
almost certainly is not. 

18026 —2415 = RST3149 AB: This pair has been seen only twice since 1935. The occultation 
pair of lAfricano et al.l (|l975l ) was not seen. 



18061-2412 = B 376: Unconfirmed since its di scovery in 1 927 (|van den Boslll928l ). The A 
component is an eclipsing binary with a 4.6 d period (|Qterd 120071 ). 

20035+3601 = STF2624Aa-B: The close pair, MCA 59Aa, which is notoriously difficult to 
detect, was not seen here. 

20181+4044 = STF2666Aa-B: The close speckle pair CHR 96Aa was not recovered. 

20205+4351 = TRN 29AC: It is unclear whether there are two new components or a nearby 
optical pair with a large relative proper motion detected at two different positions. 

21390+5729 = BU 1143AB: The close Aa pair (MIU 2) was not recovered, although it may 
have been too close. 

21449+6228 = TRN 33AC: The visual pair, MLR 16, would have been much too wide (17'^6) 
to detect here. 

22393+3903 = TRN 36AC: The visual binary, S 813, at ~ 1 arcminute of separation is 
much too wide to detect here. 

22568+6244 = TRN 37AC: It is hard to imagine a geometry where both this pair and 
MLR 266 would be dynamically stable. More than likely one (or both) are optical. 



3.1.2. Comparison ofV- and I-band Magnitude Differences 

Among the O star primaries listed in the WDS, we may estimate the companion spectral type 
from the magnitude difference l\V for those with small magnitude differences. As the magnitude 
difference gets larger, the uncertainty in the companion spectral type increases. However, the 
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situation improves with the availabihty of magnitude differences for two filters. We compare the V- 
band magnitude differences from the WDS with the I-band values we determined by AO. The result 
is shown in Figure [3j Error bars for the visual data were based on typical Tycho-2 photometric 
errors when available (smaller- Am pairs), and on scatter among individual visual estimates for those 
(larger-Am) systems with more than one such estimate. A rough average value of this scatter (0.7 
mag) was used for systems with only one visual Amy estimate. Mean Am/ values and error bars 
are given for systems with multiple observations. For systems with no /-band error estimate a 
default value of 0.5 magnitudes was adopted. 

This figure shows that AV ~ AI for small magnitude differences, but AV increases more 
rapidly at larger values of AI. This is as expected for our set of 0-star targets: the components of 
small-Am systems are of similar spectral type, but the secondaries become later in spectral type 
— hence redder — as Am increases. 

To determine the expected extent of this trend, absolute V and / magnitud es for main sequence 



spectral types were extracted from Allen's Astrophysical Quantities (|Coxll200ll ). Values of AV and 
AI were then determined, assuming primaries of spectral types 03, 05, 08, and BO and secondaries 
covering a wide range of spectral types. The results, plotted as dotted lines in Figure [3l appear 
to be in good agreement with the trend seen in these AO observations and are consistent with the 
assumption that the companions are physical. 



3.1.3. Orbits of BU 1032 (a Ori AB) and SEE 322 (HD 155889) 



Worleyl ()1990l ) provides guidelines as to when a collection of binary star measures merits the 
publication of a new orbit. Following these criteria, we determined that two systems, BU 1032 and 
SEE 322, d eserved new orbital solutions. All available measures for these pairs were extracted from 
the WDS (IMason et al.ll2001bl ) and individual measures were weighted following the precepts of 
Hartkopf et al.l ([20011). The orbital element s were determined with an iterative three-dimensional 
grid-search algorithm (jSevmour et al.l |2002| ) . For these two pairs , orbital elements are provided 
in Table [2] and future ephemerides in Table [3l Elements from the iHartkopf et al.l (119961 ) orbit for 
BU 1032 are provided for comparison. Due to the preliminary nature of the SEE 322 orbit, errors 
are not quoted. These orbits are illustrated in Figures H] and O 



Since the publication of the orbital elements of BU 1032 bv IHartkopf et al.l (|l996l ). the number 
of measures has increased by about 35%, and while the prior orbit adequately fits these data, a new 
calculation at this time makes a significant improvement in the orbit quality, as illustrated by the 
decrease in formal errors for each of the elernents. The distance and mass of the a Ori AB system 
are discussed in a recent paper by ICaballerd ([200a). 



There is no previous determination of orbital elements for SEE 322. While the errors and 
residuals are large, the orbital solution looks promising, although it is still preliminary since only 
a limited portion of the orbit is covered by existing data. The total mass associated with this 
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2 4 6 8 

Magnitude Difference (V) 

Fig. 3. — F-band magnitude differences from the WDS versus /-band values from AO. The sohd 
line indicates a 1:1 relation between AV and A J, while the various dotted curves illustrate the 
expected trend assuming main sequence primaries of types 03, 05, 08, and BO (from top to bottom, 
respectively) and a range of later spectral-type secondaries. 
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solution is very large (260 Mq for a distance of 1.2 kpc; iMason et al.lll998l ). which may indicate 
that revisions are required in the assumed elements and distance, or that the system contains more 
than two stars. 



3.2. Single Star Detections 

In addition to discovering new companions and measuring some parameters of previously 
known stars, we have a number of stars determined to be single under the specific observing 
conditions, AO system performance, and field of view. Table S] lists these stars. In this table, the 
respective columns give the HD number, Hipparcos number, epoch of observation, full-width at 
half- maximum of the corrected stellar PSF, and the V/C code for the spectroscopic binary status 
(see § 13. ip . A larger FWHM is indicative of poorer seeing, a dimmer primary, non-optimal AO 
system tuning, or some combination of these factors. 



3.2.1. Notes to Unresolved Systems 



HP 15137: This runaway star is a probable, single-lined spectroscopic binary (jBoyajian et al 
2OO5I). 



HD 25638: The pair ES 2603 (p ~ 6'.'5, Am ~ 6.0) was not detected. It was probably just out 
of the field of view. 



HD 30614: Classified as a runaway star in lGies &: BoltonI ( 19861 ) . 



HD 36879: This star is a runawav object according to its proper motion (jMdzinarishvili &: Chargeishvih 
20051 ). 



HD 37042: The known components to this st ar are all far too wide fo r detection here. The star 
9"^ Ori B is radial velocity variable according to lMorrell &: Levatd ( 199ll ). 



HD 37043: Iota Ori is a close speckle pa i r of g reat interest as, despite its very close separation, 
it was deemed optical in iGualandris et al.l (|2004l ) . The speckle pair was last measured very close 
to the FWHM value, so it is possible that the component was just barely too close to be detected. 



HD 37366: This star is a double-lined spectroscopic binary (jBovaiian et al.ll2007bl ). 



HD 39680: Classified as a runaway star in lGies &: BoltonI (|l986l ). The known pair, S 502, is 
much too wide (46'.' 1) for detection here. 

HD 41161: The known pair, ES 1234AB, was too wide (10'.'3) for detection here. Whether this 
pair is bound or not is not yet certain. 



HD 45314: CHR 251 was last measured closer (54 mihiarcsec) than the FWHM of the PSF of 



WDS 05387-0236 
P = 156.7 yr 
a = 0.2662 arcsec 



BU 1032AB 
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Fig. 4. — Relative visual orbit of BU 1032; the x and y scales are in arcseconds. The solid 
curve represents the newly determined orbit (Table [2]). The dot-dash line indicates the line of 
nodes. High-resolution measurements are shown as filled circles (speckle) or filled stars (AO). 
Visual measurements are denoted with plus signs. The ori entation and directioii of motion are 
indicated in the lower right corner of the plot. The orbit of iHartkopf et al.l (|l996l ) is shown as a 
dashed curve, which matches the new orbit quite well. 
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Fig. 5. — Relative visual orbit for SEE 322 in the same format as Fig. [4l The H indicates the 
measure of Hipparcos. 
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our ob servation. The A component is an Oe star and a probable velocity variable ([Bovaiian et al 



2007a|). 



HD 52266: This star is a probable single-lined spectroscopic binary (jMcSwain et al.l 120071 ) . 



HD 54662: This star is a double-lined spectroscopic binary (jBovaiian et al.ll2007bl ). 

HD 608 48: This Oe star show s emission line variations but no evidence of significant velocity 
variability (iBoyajian et al.ll2007al ). 



HD 69648: Not included in lMason et alj (|l998l ) survey. 

HD 163892: These data were taken using an 800-900 nm filter. 



HD 164794: This star is a probable long-period spectroscopic binary (jRauw et al.ll2002l ). 



HD 167771: The known companion, RST3170, last measured in 1940 at 8'.'3 (jRossiterl Il955l ). 
would have been outside the field of view. 

HD 175876: The known, probably optical, pair, HO 271, would have been too wide to detect. 



HD 1919 78: Not included in lMason et al.l (|l998l ) survey. The star shows no evidence of velocity 
variability (jAbt et al.lll972l ). 

HD 193443: The close pair, A 1425AB (p ~ O'.'l) may have been too close to resolve, while the 
wider AC pair (p 9'.'1) may have been outside the field of view. 



HD 195592: This is a probable single-lined spectroscopic binary ([McSwain et al.l 120071 ). 



HD 216898: Not included in iMason et al.l (|l998l ) survey. 



3.3. 2MASS Data Mining Confirmations 



Searches were made for 2MASS (jCutri et al.l 120031 ) companions to all stars on our observing 
list, using Vizier and Aladin, as well as the 'data mining' technique of lWvcoff et al.l (|2006l ). Results 
are given in Tableland discussed below. 

It should first be noted that the separation/ Am regimes covered by these two techniques have 
rather limited overlap. The "plate scale" of the AO detector nominally limits us to companions 
within ~ 5" of the primary, although this outer limit may extend to as wide as 7'.' 5, depending on 
orientation of the pair relative to the detector, as well as slight off-center placement of the primary. 

In order to get an idea of the separation/ Am limits of 2MASS, the point-source catalog was 
searched for sources in the magnitude range 5.5 < J < 8.0, corresponding to the approximate 
2MASS J-magnitude range for the AO targets in this project. This yielded 99,656 sources. All 
sources within 10" of these "primaries" were then extracted from the catalog, yielding 9,657 com- 
panions. A plot of separation versus A J is shown in Figured As seen in this figure, very few close 
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pairs are detected (only 17 pairs with separations in the range 1" < p < 2" and an additional 70 
in the range 2" < p < 3"). For the main body of companions, maximum AJ increases to perhaps 
6-6.5 mag at 3", 8-8.8 mag at 4", and 9.5-10 mag at 7'.'5. 

AO Single Stars: A Vizier search of the regions around each of the AO single stars found only 
six 2MASS companions within 7'.'5 of their respective primaries: HD 25638 (6'.'2), HD 52266 (7'.'l), 
HD 151515 (6'.'2), HD 152219 (7^2), HD 163892 (6'.'4), and HD 210809 (5'.'6). A further check of 
their orientations showed that all of these companions fell either outside the observing window or 
so close to the edge of the field as to not be measurable. The companion to HD 25638 is known; 
see S [3Xn 

Known Pairs: A search of the regions around these primaries found six 2MASS companions within 
7'.'5. Two of these (HD 48279 and 152408) correspond to known companions also measured by AO 
and are listed in Table O which lists relative astrometry and differential photometry in the 2MASS 
J, H, and Ko, bands. The remaining four, including one known companion and three objects 
much wider than the known secondaries, all fall outside the AO observing window. Note that the 
companion to HD 152408 is very red (Am/ > ArriKg), perhaps indicating that the companion still 
possesses a luminous disk (like the case of the IR-bright companion of a Ori AB). 

New Pairs: The search around these objects found five 2MASS companions within 7''5. Two of 
these (HD 156212 and 201345) correspond to the newly discovered AO companions and are listed 
in Table [5l The other three — two of which are known companions — are all wide. 



4. DISCUSSION 

Our AO survey of the 0-stars has revealed a large number of new and generally faint compan- 
ions. We show in Figure [7] the distribution of numbers of companions as a function of magnitude 
difference A/ (in 1 mag bins), and the distribution peaks near A/ = 10 mag. It is certainly pos- 
sible that some of these faint companions are low mass, gravitationally bound stars. For example, 



Hillenbrand! (|l997l ) made a deep /-band survey of the Orion cluster centered on the Orion Trapez- 
ium, and she found examples of still embedded pre-main sequence stars with magnitude differences 
as large as A/ = 15 mag compared to the bright 0-stars in the Trapezium. On the other hand, 
her survey also demonstrates that there will be many cluster stars that appear projected on the 
sky at positions near 0-stars but are not orbitally bound to them. 



We used the I-band star counts from iHillenbrandl (jl997l ) to model the possible confusion limit 



caused by nearby cluster stars that we might expect to find in our observations. The cumulative 
number distribution of stars brighter than I in her sample (1551 stars with measured I magnitudes) 
is well matched by a power law, 

N{< I) = 0.063 10^^, (1) 
where the exponent is 7 = 0.27 (valid down to / = 16 mag). The total surface density of stars in 
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Fig. 6. — Magnitude difference versus separation from 2MASS and AO. AO observations are shown 
as squares (for previously-known pairs) or diamonds (for new pairs) , with error bars plotted for the 
photometry only. The faint cloud of plus signs indicate all 2MASS pairs under 10" in separation, as 
described in the text. Two known and two new AO pairs were also found in the 2MASS data; these 
are shown as filled symbols, with 2MASS J-band magnitude differences in grey and AO magnitude 
differences in black. F-band magnitude differences for the two known pairs with AO and 2MASS 
data were very similar to the AO Am values, so are not plotted. 
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Fig. 7. — Number of companions vs. differential magnitude at /-band. Each bin represents a 
magnitude difference of 1 mag. The shaded region represents previously known companions while 
the solid line histogram shows the total number of companions (known plus new). 
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th e central part of c luster where the 0-stars reside is 10^'^ stars per square degree (see Figure 4a 



m 



Hillenbrandlll997l l or 1.9 stars per 10" x 10" field of view (= AEOS FOV). Thus, the cumulative 



distribution of cluster background stars per AEOS FOV is approximately given by 

N{< I) = 7.8 X 10"^ 10'^-^. 



(2) 



If we suppose that all 0-stars ha ve a color index of V—I = —0.39 (IWegneril 19941 ) , the extinction 
to Orion i s approximately A j = 0.25 ( Hillenbrand 1997 : Fitzpatrick 19991 ). and the distance to Orion 



is 470 pc (iHillenbr andll 1 99 71 ) . then the numbers per AEOS FOV can be recast in terms of an 0-star 
absolute magnitude My and magnitude difference A/ as 



N{< /) = 0.023 io-/i^i+Mv)^ 



(3) 



If we further assume that the same relationship holds for other 0-stars in more distant clusters, 
then the number of cluster stars appearing in the FOV will scale as distance d squared. Let r/ 
represent the predicted number of cluster background stars per AEOS FOV. The relation between 
the limiting magnitude difference A/ associated with ry is then 



2 1 
A/ = 3.6 log d - My - log rj. 

7 7 



(4) 



Thus, we expect that background confusion with cluster stars will increase (smaller A/) with 
increasing distance and with lower luminosity 0-stars. 

We applied this relation to estimate a limiting A/ for each cluster and association star in our 
sample in order to separate the probable physical comp anions from the fai nt optical companions 
caused by background confusion. We used distances f rom iMason et al.l (jl998l ) and the spectral type 
- absolute magnitude calibration from lMartins et al.l (j2005l ). The Orion cluster richness, extinction, 
and age may not be representative of the full sample of cluster 0-stars in our survey, and so the 
normalization of the A/ relation is somewhat uncertain. We note that large numbers of companions 
begin to appear near AI = 8 mag (see Figure [7]), and for typical 0-star distances (d = 1 to 2 kpc) 
and magnitudes (My = —4 to —6) this corresponds to a background density in the Orion model of 
r] = 1 star per AEOS FOV. Thus, we adopted this value of rj in equation H] to identify the probable 
physical companions among the cluster and association stars. 

The background limits for the field and runaway stars depend instead upon on the stellar 
number density of Galactic background stars as a function of magnitude and the direction of a 
given target. Here we adopted the cumulative star counts as a fui iction of apparent / magnitude 
and Galactic coordinates from the Besangon model of the Galaxy (jRobin et al.l 120031 ). We created 
a cumulative star counts versus AI relation for each of these stars and found a limiting magnitude 
for a specified value of 0.01 stars per AEOS FOV. We set this conservative limit based on numerical 
experiments for the set of runaway stars with observed faint companions. None of these companions 
can be physical companions since any such wide binary would be disrupted at the time the runaway 
star was ejected (by close gravitational encounters or a supernova in a binary). 
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The results of this exercise to discern the probable physical companions are summarized in 
Table [6j The targets were placed in categories of cluster and associa tion stars, field stars, and 
runaway stars according to the criteria adopted bvlMason et alj (|l998l ) (with some revisions from 
the subsequent work on field 0-stars bv lde Wit et al.ll2004l . l2005l ). We see that after rejection of the 
possible background stars there is still a high frequency (37%) of detected companions among the 
cluster and association stars in our survey. In sharp contrast, there are very few such companions 
amo ng the field a n d ruii away stars. The only field star with a probable companion is HD 48279, 
and Ide Wit et al.l (|2004l ) show that this star is 28" from a cluster identified in 2MASS images, 
which suggests that this 0-star may be a member of a young, emerging cluster. The only runaway 
star with a candidate co mpanion is HD 34078 = AE Aur, a star ejected from the Orion association 



Gualandris et al 



Mason et al 



20041 ). There is another optical companion, SEI 136, at a separation of 8'.'4 



19981 ). and we suspect that both the companion we found and SEI 136 form a visual 
system that appears by a chance alignment with AE Aur at this time in its trajectory across the 
sky. 



Our results confirm the trends found by lMason et al.l (119981 ) that binary 0-stars are common 
among 0-stars found near their birthplaces and are rare among 0-stars ejected from clusters. These 
trends are consistent with curre nt models for the ejection of runaways and the formation of massive 
stars (jZinnecker &: Yorkd 120071 ). High velocity 0-stars were probably ejected b y close encounters 



betw een binaries in dense clusters and by supernova explosions in close binaries (jHoogerwerf et al 



2001 



), and their ejection velocities generally exceed the escape velocity binding any wide, orbiting 
companion. Many massive stars are probably born ii i high density ste llar environments where bina- 



ries may be formed through three-body i nteractions (IBate et al 



with large proto-stellar accretion disks (iMoeckel &: Ballyi 120071 ) . Even in lower de nsity environ 



2002) and through the interactions 



ments, binaries may form through disk fragmentation and subsequent gas accretion (iKratter et al 



20071 ). Our results add to growing body of evidence that the formation of binaries is closely linked 
to the formation of massive stars and the deposition of the angular momentum reservoir of the 
natal cloud. 
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Table 1. Survey Measurements 



Discoverer Epoch P.A. Sep. Am/ Spec. 



WDS 


HD 


HIP 


Designation 


(BY) 


(°) 


(") 




(mae) 




Stati 


(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


(7) 




(8) 




(9) 


00061+6341 


108 


505 


TRN 


7 




2001, 


.7450 


5Z/176N* 


3.24 


9, 


.47+0, 


.26 


C 


00177+5126 


1337 


1415 


TRN 


8 




2001, 


.7341 


22Z/186N'^ 


2.75 


9, 


,30+0, 


,47 


V 


02158+5600 


13745 


10541 


TRN 


9 




2004, 


,7924 


29 


4.55 


10, 


.04+0, 


.67 


V 


02229+4129 


14633 


11099 


TRN 


10 




2001, 


.7451 


97Z/214N^ 


2.18 


10, 


.31+0, 


.61 


V 


02327+6127 


15558 


11832 


TRN 


11 


AF 


2001, 


.7450 


261 


4.45 


6, 


.59+0, 


.19 


V 














2002, 


.6845 


261 


4.48 


6, 


,38+0, 


,40 
















2005, 


.6962 


262 


4.44 


6, 


.42+0, 


.30 










TRN 


11 


AG 


2001, 


.7450 


221 


2.99 


9, 


.94+0, 


.35 










TRN 


11 


AH 


2001, 


.7450 


211 


5.79 


11, 


.03+0, 


.86 




02407+6117 


16429 


12495 


TRN 


12 


AD 


2005, 


.6963 


112 


2.94 


7, 


,39+0, 


,22 


V 


02511+6025 


17505 


13296 


STF 


306 


AB 


2001, 


.7450 


91 


2.09 


1, 


.71+0, 


.26 


V 














2002, 


.6845 


89 


2.11 


1, 


.71+0, 


.37 
















2004, 


.9648 


92 


2.10 


1, 


,64+0, 


,37 










TRN 


13 


AH 


2001, 


.7450 


142 


4.59 


6, 


,72+0, 


.21 
















2002, 


.6845 


140 


4.62 


7, 


.01+0, 


.27 
















2004, 


.9648 


143 


4.59 


7, 


.15+0, 


.22 




02594+6034 


18326 


13924 


TRN 


14 




2001, 


.7451 


223Z/353N^ 


2.38 


7, 


,45+0, 


,46 


V 


03141+5934 


19820 


15063 


TRN 


15 




2001, 


.7452 


81Z/209N'^ 


2.81 


10, 


.37+0, 


.42 


V 


03556+5238 


24431 


18370 


HDS 


494 




2001, 


.7455 


176 


0.71 


2, 


.86+0. 


.35 


c 


03590+3548 


24912 


18614 


TRN 


16 




2001, 


.7342 


141Z/236N^ 


2.40 


9, 


.81+1. 


.16 


c 


05163+3419 


34078 


24575 


TRN 


17 


Aa 


2004, 


.7873 


171 


0.35 


4, 


.97+0, 


,51 


c 


05207+3726 


34656 


24957 


TRN 


18 


Aa 


2002, 


.7257 


279 


0.35 


3, 


.77+0, 


,33 


c 














2003, 


.7360 


280 


0.35 


4, 


.09+1, 


,10 
















2003, 


.8917 


280 


0.34 


4, 


,13+0, 


,84 










TRN 


18 


AB 


2002, 


.7257 


48 


1.88 


6, 


,98+0, 


,24 
















2003, 


.7360 


47 


1.90 


6, 


.94+0, 


,25 
















2003, 


.8917 


47 


1.89 


7, 


.28+0, 


.43 










TRN 


18 


AC 


2003, 


.7360 


126 


2.96 


10, 


.57+0, 


.68 










TRN 


18 


AD 


2003, 


.7360 


247 


5.75 


10, 


.39+0, 


.42 




05297+3523 


35921 


25733 


HU 


217 




2004, 


.7873 


253 


0.61 


2, 


.23+0. 


.31 


V 


05351+0956 


36861 




STF 


738 


AB 


2001, 


.7428 


297Z/356N'*'' 


'4.25 


2, 


.29+0. 


.16 


c 














2001, 


.7456 


41 


4.23 


2, 


,04+0, 


,15 




05353-0523 


37022 


26221 


STF 


748 


Ca,F 


2001, 


.7456 


118 


4.39 


4, 


.89+1, 


,06 


V 














2003, 


.0045 


120 


4.46 


4, 


.73+0. 


.19 




05354-0525 


37041 


26235 


CHR 


249 


Aa 


2004, 


.8666 


292 


0.39 


2, 


.97+0. 


.72 


V 
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Table 1 — Continued 









Discoverer 


Epocli 


P.A. 


Sep. 


Ami 


Spec. 


WDS 


HD 


HIP 


Designation 


(BY) 


n 


(") 


(mag) 


Status'^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



05387-0236 


37468 


26549 


BU 


1032 AB 


2001, 


.7456 


106 


0, 


.24 


1. 


.22+2, 


.15 


V 








TRN 


19 AF 


2001, 


.7456 


18 


3 


.13 


8, 


,07+0, 


,33 




05407-0157 


37742 




STF 


774 Aa-B 


2001, 


.7455 


160 


2, 


.36 


2, 


.24+0, 


.89 


C 


06319+0457 


46150 


31130 


GAN 


3 AB 


2001, 


.0989 


284 


3 


.47 


4, 


.87+0, 


.13 


C 


06322+0450 


46223 


31149 


TRN 


20 


2001, 


.8661 


314Z/264N* 


0, 


.46 


4, 


.37+0, 


,17 


c 


06364+0605 


46966 


31567 


TRN 


21 


2002, 


.8023 


210 


3 


.19 


9, 


,97+0, 


,63 


c 


06374+0608 


47129 


31646 


TRN 


22 AB 


2002, 


.2404 


240 


1, 


.15 


5, 


,10+0, 


,43 


V 












2005, 


.0144 


251 


1, 


.16 


4, 


.94+0, 


,22 










TRN 


22 AC 


2005, 


.0144 


203 


0, 


.78 


5, 


,14+0, 


,22 




06386+0137 


47432 


31766 


TRN 


23 


2005, 


.0134 


206 


0, 


.78 


4, 


,95+1, 


,31 


c 


06410+0954 


47839 


31978 


STF 


950 Aa-B 


2001, 


.7457 


202 


2, 


.88 


3, 


.31+2, 


,11 


V 


06427+0143 


48279 


32137 


STF 


956 AB 


2004, 


.0478 


199 


6, 


.63 


2, 


.57+0. 


.39 


c 


08392-4025 


73882 


42433 


B 


1623 


2004, 


.1190 


256 


0, 


.65 


1, 


,25+0, 


,77 


V 


16550-4109 


152408 


82775 


I 


576 


2001, 


.4030 


260Z/261N'' 


5, 


.30 


5, 


.67+0, 


,92 


c 


17065-3527 


154368 


83706 


B 


894 


2002, 


.6730 


357 


2, 


.57 


6, 


.28+0. 


.54 


V 


17158-3344 


155889 


84444 


SEE 


322 


2003, 


.6041 


285 


0, 


.19 


0, 


.68+2. 


.59 


c 


17175-2746 


156212 


84588 


TRN 


24 AB 


2001, 


.4959 


206Z/204N'^ 


4, 


.04 


8, 


,05+0, 


,39 


c 








TRN 


24 AC 


2001, 


.4959 


145Z/144N'^ 


7, 


.31 


7, 


.40+0, 


,24 




17347-3235 


159176 


86011 


HDS 


2480 Ab 


2003, 


.3280 


59 


0, 


.70 


3, 


.14+0, 


,47 


V 


17595-3601 


163758 




TRN 


25 


2001, 


.4960 


107Z/ 93N^ 


1, 


.70 


8, 


,05+0, 


,33 


c 


18024-2302 


164492 


88333 


TRN 


26 AH 


2002, 


.4516 


342 


1, 


.48 


5, 


,12+0, 


,23 


V 


18026-2415 


164536 




RST 


3149 AB 


2002, 


.4627 


62 


1, 


.65 


4.49+12'i 


c 


18061-2412 


165246 




B 


376 


2002, 


.4819 


98 


1, 


.90 


3, 


,68+0, 


.47 


V 


18152-2023 


167263 




BU 


286 AB 


2005, 


.6545 


216 


5, 


.93 


5, 


,62+0, 


,98 


V 


18181-1215 


167971 


89681 


TRN 


27 


2001, 


.7447 


80Z/ 39N'^ 


4, 


.70 


8, 


.12+0, 


.32 


V 


20035+3601 


190429 


98753 


STF 


2624 Aa-B 


2001, 


.7337 


172 


1, 


.92 


0, 


.84+0. 


.95 


c 


20181+4044 


193322 


100069 


STF 


2666 Aa-B 


2001, 


.6683 


244 


2, 


.69 


2, 


.51+0. 


.12 


V 












2001, 


.7364 


243 


2, 


.67 


2, 


,11+0, 


,26 




20191+3916 


193514 


100173 


TRN 


28 AB 


2001, 


.7367 


78 


4, 


.75 


5, 


,80+0, 


,17 


c 












2002, 


.6730 


79 


4, 


.79 


6, 


.18+0. 


.25 










TRN 


28 AC 


2001, 


.7367 


157 


3, 


.20 


7, 


,23+0, 


,28 














2002, 


.6730 


158 


3 


.22 


7, 


,44+0, 


,16 




20205+4351 


193793 


100287 


BU 


1207 AB 


2001, 


.3400 


209 


4, 


.77 


8, 


,44+0, 


,31 


V 












2001, 


.7337 


209 


4, 


.77 


6, 


.75+0. 


.28 














2001, 


.7363 


209 


4, 


.77 


6.75+12'^ 
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Table 1 — Continued 









Discoverer 


Epoch 


P.A. 


Sep. 


Am/ 


Spec. 


WDS 


HD 


HIP 


Designation 


(BY) 


n 


(") 


(mag) 


Status^^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 













2001.7447 


b 


4.79 


6.43+0.53 














2002.6731 


210 


4.81 


6.58+0.13 














2002.6731 


210 


4.81 


6.93+0.12 










TRN 


29 AC 


2001.7363 


202 


2.29 


9.62+0.35 














2001.7447 


b 


3.23 


9.32+0.38 




20566+4455 


199579 


103371 


TRN 


30 


2001.7338 


30Z/205N^ 


3.76 


9.79+0.29 


V 


21079+3324 


201345 


104316 


TRN 


31 


2001.7367 


342Z/218N'^ 


7.38 


9.31+0.47 


C 


21185+4357 


203064 


105186 


TRN 


32 


2001.7337 


81Z/219N'^ 


3.84 


9.48+0.29 


c 


21390+5729 


206267 


106886 


BU 


1143 AB 


2001.7338 


317 


1.78 


5.65+0.21 


V 


21449+6228 


207198 


107374 


TRN 


33 AC 


2001.7364 


75Z/207N'* 


2.96 


10.43+0.52 


c 


22021+5800 


209481 


108772 


TRN 


34 


2001.7341 


358Z/208N'^ 


2.75 


9.92+0.46 


V 


22051+6217 


209975 


109017 


TRN 


35 AD 


2001.7367 


341Z/171N^ 


4.14 


9.98+0.44 


c 








TRN 


35 AE 


2001.7367 


344Z/175N'* 


3.79 


10.03+0.37 




22393+3903 


214680 


111841 


TRN 


36 AC 


2001.7339 


85Z/222N'^ 


3.56 


9.94+0.45 


c 


22568+6244 


217086 


113306 


MLR 


266 AB 


2001.7369 


353 


2.79 


3.41+0.27 


c 












2002.6843 


354 


2.83 


3.34+0.48 














2004.7598 


354 


2.78 


3.79+0.58 










TRN 


37 AC 


2001.7369 


164 


3.15 


7.14+0.29 














2002.6843 


165 


3.17 


6.79+12'i 














2004.7598 


165 


3.12 


7.03+0.34 





^Image dcrotator in an unknown state: Z = zenith up and N = north up. Confirmation of the pair will 
clear up this residual ambiguity. 

''Image derotator in an unknown state. Neither Zenith up or North up match coordinates of known object 
in field. 

'^Codes are: V - one or more additional spectroscopic companions indicated, or C - radial velocity constant. 
■^This is a lower limit for the errors, the actual errors are typically much larger. 



Table 2. Orbital Elements 



WDS Desig. 


Discoverer 


Period 


Semimajor 


Inclination 


Longitude 


Epoch of 


Eccentricity 


Longitude of 




Designation 




Axis 




ot IN ode 


Periastron 




Periastron 


a5 (2000) 




P (yr) 


«(") 


i(o) 


n (o) 


To (yr) 


e 


UJ (o) 


05387-0236 


BU 1032 AB 


156.7 


0.2662 


159.7 


121.7 


1999.5 


0.0515 


8.7 






±3.0 


±0.0021 


±3.7 


±9.6 


±10.2 


±0.0080 


±16.9 


05387-0236 


BU 1032 AB^ 


155.3 


0.2642 


160.4 


136. 


1997. 


0.051 


18. 






±7.5 


±0.0052 


±7.2 


±25. 


±24. 


±0.015 


±37. 


17158-3344 


SEE 322 


536. 


0.35 


111. 


167. 


1870. 


0.64 


109. 



'^Elements from 



Hartkopf et al. 
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Table 3. Future Ephemerides (BY) for BU 1032 and SEE 322 



WDS Desig. 


Discoverer 


2008 


2010 


2012 


2014 


2016 


a, 5 (2000) 


Designation 


e n p (") 


e n p (") 


e n p (") 


6 n p (") 


e n p (") 


05387-0236 


BU 1032 AB 


93.0 0.249 


88.0 0.249 


83.1 0.248 


78.1 0.247 


73.1 0.246 


17158-3344 


SEE 322 


281.7 0.184 


280.4 0.184 


279.0 0.183 


277.7 0.183 


276.3 0.183 
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Table 4. Single Star FWHM measures 



Epoch FWHM Spec. 
HD HIP (BY) (") Status^ 



14947 


11394 


2001.7455 


0.12 


C 


15137 


11473 


2004.7923 


0.18 


V 


25638 


19272 


2004.8665 


0.12 


c 


30614 


22783 


2004.8665 


0.12 


c 


36879 


26272 


2003.8944 


0.11 


c 


37042 




2004.8667 


0.12 


V 


37043 


26241 


2001.7455 


0.08 


V 


37366 


26611 


2005.0200 


0.34 


V 


39680 


27941 


2002.2404 


0.22 


c 


41161 


28881 


2001.8659 


0.16 


c 


42088 


29216 


2001.8660 


0.21 


c 


45314 


30722 


2001.0989 


0.11 


V 


46149 


31128 


2001.8660 


0.36 


V 


46966 


31567 


2002.2377 


0.25 


c 


48099 


32067 


2002.2405 


0.13 


V 


52266 


33723 


2001.9754 


0.12 


V 


53975 


34297 


2005.0196 


0.12 


V 


54662 


34536 


2002.0218 


0.16 


V 






2002.2404 


0.12 




55879 


34999 


2005.0170 


0.19 


c 


57061 


35415 


2004.0807 


0.21 


V 


57682 


35707 


2001.9290 


0.45 


c 


60848 


37074 


2001.0990 


0.16 


c 


66811 


39429 


2004.0861 


0.48 


c 


69648 




2004.0451 


0.26 




75211 


43125 


2004.1432 


0.23 




75222 


43158 


2005.0197 


0.49 


c 


93521 


52849 


2001.0992 


0.19 


c 


148546 


80829 


2005.3482 


0.28 


c 


149404 


81305 


2001.3893 


0.12 


V 


149757 


81377 


2002.2243 


0.22 


c 
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Table 4 — Continued 



Epoch FWHM Spec. 
HD HIP (BY) (") Status^^ 



151003 


82121 


2004, 


.2590 


0. 


,23 


V 


151515 


82366 


2004, 


.3384 


1. 


,93 


c 


152003 




2001, 


.3893 


0. 


,12 


c 


152219 




2004, 


.5265 


0. 


,22 


V 


152314 




2001, 


.4029 


0. 


.17 


c 


153919 


83499 


2003, 


r" Ci Ai^ 

.5247 


0. 


,16 


V 


155806 


84401 


2002, 


.6730 


0. 


1 O 

,12 


c 


157857 


85331 


2002, 


.2383 


0. 


.22 


C 


162978 


87706 


2001, 


.7446 


0. 


,12 


c 


163758 




2001, 


.5179 


0. 


,17 


c 


163800 


88040 


2001, 


.4960 


0. 


,13 








2001, 


.5176 


o 

0. 


lb 




163892 


88085 


2005, 


.3404 


0. 


,21 


V 


1 A /loo 

164438 


88297 


2002, 


o on A 

.3204 


0. 


lb 




164794 


88469 


2002, 


.6729 


0. 


,10 


V 


164816 




2001, 


.6709 


0. 


1 A 

,14 


V 






2002, 


.4817 


0. 


,11 




165052 


88581 


2001, 


.6711 


0. 


,14 


V 


165319 


88652 


2001, 


.6711 


0. 


16 


c 


167771 


89630 


2002, 


.4843 


0. 


,20 


V 






2002, 


.4950 


0. 


,19 




175876 


93118 


2001, 


.7447 


0. 


,14 


c 


186980 


97280 


2001, 


.7366 


0. 


,19 


c 


188001 


97796 


2001, 


.7338 


0. 


,12 


V 






2001, 


.7447 


0. 


,17 




188209 


97757 


2004, 


.7646 


0. 


,08 


c 


190864 


98976 


2001, 


.7365 


0. 


,14 


V 


190918 


99002 


2001, 


.7421 


0. 


,18 


V 


191612 


99308 


2001, 


.7365 


0. 


,13 


c 


191978 


99439 


2002, 


.6786 


0. 


33 


c 


192281 


99580 


2001, 


.7365 


0. 


,12 


c 
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Table 4 — Continued 







Epoch 


FWHM 


Spec. 


HD 


HIP 


(BY) 


(") 


Status*^ 


192639 


99768 


2001.7366 


0.17 


C 


193443 


100146 


2001.7366 


0.19 


V 


195592 


101186 


2005.6549 


0.19 


V 


198846 


102999 


2005.3843 


0.12 


V 


202124 


104695 


2005.6549 


0.19 


c 


203064 


105186 


2001.7342 


0.12 


c 






2001.7448 


0.12 




207538 


107598 


2005.7180 


0.22 


c 


210809 


109562 


2001.7369 


0.21 


c 


210839 


109556 


2001.7342 


0.16 


c 


214680 


111841 


2001.7369 


0.19 


c 


216898 


113218 


2001.7451 


0.10 




218915 


114482 


2004.7597 


0.18 


c 



'^Codes are: V - one or more additional spectro- 
scopic companions indicated, or C - radial velocity 
constant. 



Table 5. Data Mining Results from 2MASS 



Discoverer 


HD 


HIP 


WDS 


Epoch 


P.A. 


Sep. 


Amj 




Am^^ 


Designation 








(BY) 


n 


(") 


(mag) 


(mag) 


(mag) 


STF 956 AB 


48279 


32137 


06427+0143 


1999.87 


194.8 


6.55 


1.70 


1.63 


1.56 
















±0.06 


±0.03 


±0.04 


I 576 


152408 


82775 


16550-4109 


1999.36 


264.2 


4.72 


1.47 


1.71 


1.25 
















±0.35 


±0.58 


±0.44 


TRN 24 AC 


156212 


84588 


17175-2746 


1998.53 


147.6 


7.32 


4.06 


5.38 


4.24 


















±0.14 




TRN 31 


201345 


104316 


21079+3324 


1999.75 


220.8 


7.44 


6.38 


5.88 


5.69 
















±0.89 


±0.45 


±0.24 
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Table 6. Companion Frequency and Environment 
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Number 


Number 


Number 


Percentage 




of Stars 


w/ Comp. 


w/ Phys. Comp. 


of Phys. Comp. 


Cluster / Association 


83 


41 
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Field 


9 


2 
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11 ± 11 
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24 


8 


1 
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